Measuring the Cosmic Equation of State with Counts of Galaxies 
II: Error Budget for the DEEP2 Redshift Survey 
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ABSTRACT 

In a previous paper, we described a new variant on the classical "dA'^/dz" test which could 
be performed using data from the next generation of redshift surveys. By studying the apparent 
abundance of galaxies as a function of their circular velocity or velocity dispersion, rather than 
luminosity, it is possible to avoid many of the uncertainties of galaxy evolution while using 
quantities that may be measured directly. In that work, we assumed that counting statistics 
would dominate the resulting errors. Here, we present the results of including cosmic variance 
and determine the impact of systematic effects on attempts to perform the test with the upcoming 
DEEP2 Redshift Survey. For the DEEP2 survey geometry, cosmic variance yields errors roughly 
twice those predicted from Poisson statistics. Through Monte Carlo simulations we find that 
if the functional form, but not the strength, of any of the major systematic effects (baryonic 
infall, velocity errors, and incompleteness) is known, the free parameter may be determined 
from the observed velocity function. The systematic may then be corrected for, leaving a much 
smaller residual error. The total uncertainty from systematics is comparable to that from cosmic 
variance, but correlated amongst redshift bins. Based on these analyses, we present error budgets 
for a dN /dz measurement with DEEP2 and determine the resulting constraints on cosmological 
parameters. We find that the uncertainty in the cosmic equation of state parameter w are ~ 2 x 
higher than previously derived, providing a measurement much stronger than any available today 
but weaker than some other proposed tests. 



Subject headings: cosmological parameters, cosmology: 
fundamental parameters 

1. Introduction 

In a previous paper (Newman & Davis 2000, 
hereafter NDOO), we described a new variant on 
the classical "dA^/dz" test which could be per- 
formed using data from the next generation of 
redshift surveys. This classical technique deter- 
mines the evolution of the cosmic volume element 
by measuring the redshift distribution of a tracer 
whose number density is known, providing con- 
straints on fundamental cosmological parameters. 
In the past, the abundance of galaxies was used to 
perform this test under the assumption that the 
total comoving number density of galaxies inte- 
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grated over all luminosities is independent of red- 
shift (with further assumptions about the lumi- 
nosity function; e.g. Loh & Spillar 1986). Such 
assumptions may be suspect, but improving upon 
them would require a reliable, comprehensive the- 
ory of galaxy formation and evolution. 

Instead, NDOO suggests that by measuring the 
apparent abundance of galaxies as a function of 
their linewidth or velocity dispersion (normalized 
to the abundance at z 0), it is possible to ob- 
tain a more reliable measurement of the volume 
element, exploiting the simplicity of the dark mat- 
ter halo velocity function. Current and upcom- 
ing large redshift surveys with moderately high 
spectroscopic resolution will make this possible. 
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The DEIMOS/DEEP Rcdshift Survey (hereafter, 
"DEEP 2") is particularly well suited for the tech- 
nique. This project is a major effort to observe 
distant galaxies using the new DEIMOS spectro- 
graph on the Keck telescope, scheduled to be in- 
stalled in early 2002. Its goal is to obtain high 
quality spectra of 50,000 galaxies selected to have 
minimum redshift z > 0.7 (the "IHS", or 1-hour 
survey, so named because of the expected expo- 
sure time per slitmask) and spectra of 5,000 such 
galaxies in selected regions to a fainter limiting 
magnitude (the "3HS", or 3-hour survey). DEEP2 
will obtain data characterizing galaxies and large- 
scale structure at z 1 comparable in quality to 
what is currently available for ^; = 0. Although 
we have much information on the local universe, 
we know little about how it has reached its present 
state; this survey has been designed to address this 
major gap in our knowledge. 

In our previous work, we assumed that the un- 
certainties in a dN/dz measurement using DEEP2 
would be dominated by Poisson statistics and 
found that it could provide strong constraints 
on the cosmic equation of state parameter of 
quintessence-like models, w — P/p. Because the 
velocity function of dark matter halos is so simple 
(a nearly perfect power-law at galactic scales), we 
presumed that any systematic effects would leave 
a clear signature and could be removed or avoided, 
and that counting statistics would therefore dom- 
inate the uncertainties. However, we did not con- 
sider the fundamental limit to the accuracy of any 
count imposed by cosmic variance. In this paper, 
we calculate its impact on the DEEP2 survey. 

Since NDOO, a mmiber of papers suggesting that 
systematic effects would compromise our proposed 
test have appeared in the literature (e.g. Bullock 
et al. 2000), though again without actually evalu- 
ating their impact on realistic methods of measur- 
ing dN/dz. This paper addresses that gap in our 
previous work; we have performed quantitative 
tests of the impact systematic effects should have 
on our ability to measure dN/dz with DEEP2. 
However, we must do so without the information 
that will be available at the completion of the sur- 
vey. For instance, semi-analytic and N-body esti- 
mates of the impact of baryonic infall on the po- 
tential well depths of galaxies are still developing 
(e.g. Kochanck & White 2001), and comparisons 
between HI rotation curves and 0[II] linewidths 



for DEEP2 galaxies may be possible in the future, 
but are unavailable to us now. 

Rather than attempting to foresee all that the 
next few years will bring, we have instead cho- 
sen to evaluate the impact of major systematic ef- 
fects by applying "toy" models of their impact to 
a semi-analytic prediction of the dark-matter halo 
velocity function (a la NDOO). The semi-analytic 
methods used duplicate the general behavior seen 
in N-body simulations. The toy models used each 
have one free parameter that is treated as un- 
known going into the analysis. We then generate 
Monte Carlo realizations of the velocity function 
affected by these model systematics and determine 
the resulting errors in measuring both the free pa- 
rameter and the abundance of dark halos in such 
datasets. Thus, for any value of the free parameter 
we can determine the excess variance in a dN/dz 
measurement due to the residual uncertainty from 
that systematic effect after removal. So long as 
the toy models duplicate the major features of the 
real systematics, we should obtain reasonable es- 
timates of the actual errors from DEEP2. In § 2 
of this paper, we calculate the impact of cosmic 
variance; in § 3, we determine the residual errors 
that will result from systematic effects; and in § 
4, we present the constraints on cosmological pa- 
rameters that a measurement of dA'^/d2; from the 
velocity function of DEEP2 galaxies would pro- 
vide. 

2. Cosmic Variance 

By definition, the spatial correlations of galax- 
ies imply that they are not distributed indepen- 
dently and thus cannot obey simple Poisson statis- 
tics. Instead, in a given volume there will be an 
excess variance 

^cv = ^j ^{xi,X2)cfxid^X2, (1) 

where ^ is the two-point correlation function of 
galaxies and V is the volume in which we wish 
to calculate the variance. This integral may more 
easily be evaluated in Fourier space, where we ob- 
tain 

4y = ^ / Pik)\W{k)\'d'k, (2) 

where W{k) is the Fourier transform of a spatial 
window function that is 1 inside the volume and 
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zero elsewhere and P{k) is the power spectrum of 
the galaxy distribution. To evaluate the integral, 
we use the CDM-like power spectrum of Bardeen 
et al. (1986) with T = 0.25 and slope n = 1, 
consistent with recent observations. 

We are interested specifically in the effect of 
cosmic variance upon measurements of diV/d2: 
with DEEP2 (as the volume observed in the 
next generation of 2; ~ surveys is much 
greater, cosmic variance in local comparison sam- 
ples should be comparatively negligible). Be- 
cause the CFH12k camera used in imaging has 
a 30' X 40' field of view, the most efficient ge- 
ometry for each of the four fields to be surveyed 
(Davis & Faber 1998) is one which is 30' in the 
shortest direction. Based on observing time con- 
straints, we further expect each field to cover 120' 
in the long direction on the sky; the required imag- 
ing is nearly complete. If we assume an rim— 0.3 
ACDM model, this corresponds to 20.2 x 80.8h~^ 
Mpc comoving a,t z = 1, the expected median 
redshift for the survey. We have calculated the 
cosmic variance for each field assuming that the 
volume surveyed is a rectangular solid extending 
1300 Mpc comoving in the redshift direction 
(corresponding to 0.7 < z < 1.5). 

For a rough estimate, we could assume that the 
power spectrum of galaxies has the same normal- 
ization as clusters today, applying the as measure- 
ment of Borgani et al. (1999), which yields 0.96 
for an flm — 0.3 LCDM model (for quintessence 
models, we use the relations of Wang and Stein- 
hardt (1998) with 6 = 0.075). We then find that 
the uncertainty in a density measurement will be 
5.9% in each of the four fields, or 2.9% for the en- 
tire survey (as the fields are widely spaced on the 
sky, there should be minimal correlation between 
them at 2: = 1, so the combined variance goes as 

1/ N fields)- 

However, to calculate the variance in counts of 
galaxies, we must use the power spectrum of galax- 
ies &t z = 1 rather than clusters at 2; = 0; we 
therefore must renormalize to match cts for galax- 
ies. A measurement of the correlation length ro 
may be directly transformed into a measurement 
of (Tg if the correlation index 7 is known (Pea- 
cock 1999); however, the correlation properties of 
galaxies at 2 ~ 1 remain poorly known. Obser- 
vations suggest that ro ~ 2.5h~l Mpc (comov- 
ing; see Hogg et al. 2000), similar to the cor- 



relation length of the dark matter at that red- 
shift, but semi-analytic models predict ro ~ 4/i~^ 
Mpc (comoving) for luminous galaxies at that 
redshift (Benson et al. 2001). For 7 = 1.8 and 
ro = 2.5,3.25,4, as = 0.474,0.600,0.723. Thus, 
we expect cosmic variance in the DEEP volume 
to be somewhat lower than a naive calculation 
would yield: 1.4 - 2.2 % rather than 2.9%. Since 
DEEP2 will obtain linewidths for ~ 10^ galax- 
ies with linewidths, the uncertainty due to cosmic 
variance will be roughly twice the Poisson value. 
In Figs. 1 and 2, we show the effect of changing the 
survey geometry or arcal coverage on the cosmic 
variance in overall DEEP2 results. The amount 
of area surveyed on the sky is clearly of greater 
importance than the shape of the surveyed region. 

Cosmic variance need not impose a permanent 
limit to the precision of a DEEP2 dN/dz mea- 
surement; if imaging of a wider area is available 
sufficient to produce an accurate photometric red- 
shift distribution (even if information on individ- 
ual galaxies remains uncertain), it could be possi- 
ble to renormalize the total abundances of galaxies 
in the DEEP2 fields and reweight the abundances 
of galaxies with observed linewidths. Other red- 
shift surveys extending to 2 ~ 1, particularly the 
VLT/VIRMOS survey, could also provide correc- 
tions. However, cosmic variance would almost cer- 
tainly remain an insurmountable obstacle to ef- 
forts to measure the evolution in the equation of 
state parameter w via DEEP2 dA''/d2 alone, since 
the variance is higher in any subvolume than for 
the survey as a whole. 

3. Correction for systematic effects 

Like most astronomical studies, the measure- 
ments proposed by NDOO could be subject to a 
number of systematic effects. If these systematics 
are well-understood and their magnitude is known, 
it is straightforward to correct for them when mea- 
suring the halo velocity function. However, few, if 
any, of the systematic effects that will affect this 
measurement can yet be predicted accurately. For 
instance, we can determine a priori what effect 
a given luminosity-linewidth relation would have 
upon incompleteness in the velocity function, but 
without sufficient measurements at 2 ^ 1 (which 
will not be possible until the DEEP2 survey is con- 
ducted), we cannot determine the resulting uncer- 
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tainties with precision. 

Thus, rather than attempt to guess the magni- 
tude of systematic effects that will be better un- 
derstood five years from now and calculate the er- 
rors in a measurement of dA''/d^; that would result, 
we have attempted to determine to what degree we 
can measure and correct for systematics using the 
signatures they leave on the velocity function at 
0^1. In many cases, this is a pessimistic as- 
sumption; for instance, we will be able to measure 
our errors in determining linewidths by compar- 
ison to absorption line measurements of DEEP2 
galaxies, rather than having to infer them from 
the velocity function alone. However, this tech- 
nique has the significant advantage that we do not 
need to understand perfectly either the intrinsic 
velocity function of dark halos or any of the sys- 
tematic eff'ects to be able to assess their impact, so 
long as we can define "toy" models for them that 
replicate their major features and have the same 
free parameters that must be determined from the 
observations. For the intrinsic dark matter veloc- 
ity function, we use the semi-analytic method of 
NDOO. We have tested models of three major sys- 
tematic effects expected in our measurements: the 
impact of baryonic infall on the circular velocities 
of dark halos; the effects of incompleteness at low 
circular velocities; and random errors in determin- 
ing those circular velocities from observations. We 
describe our models for these in more detail below. 

As a further simplification, since all physical 
models of either the halo velocity function or sys- 
tematic effects simultaneously predict their depen- 
dence upon redshift, we determine the impact of 
systematics on AN /Az measurements by treating 
DEEP2 as simply measuring the velocity function 
at z — 1 (the expected median redshift of the sur- 
vey). This should provide the same constraints 
that would result from likelihood maximization 
over redshift and velocity simultaneously, but re- 
duces the number of variables to be evaluated and 
thereby speeds computation. Thus, to determine 
the uncertainty resulting from one of the three sys- 
tematic effects we: 1) Choose a "true" value for 
the free parameter a of the systematic, while us- 
ing fiducial values for the other two. This defines a 
single probability distribution function (PDF) for 
galaxies (which for this purpose we define as all 
objects with Vc < 300 km s~^). 2) Generate 500 
Monte Carlo realizations of the velocity distribu- 



tion of DEEP2 galaxies by drawing an appropri- 
ate number of objects (e.g. 10,000 for our "best 
bet" scenario) from this PDF. 3) For each real- 
ization, we attempt to determine the input value 
of the free parameter by likelihood maximization 
(using PDFs defined on a dense grid in a). The 
distribution of the best-fit values from these real- 
izations defines the error in measuring a; this may 
be propagated into the uncertainty in the dN/dz 
measurement. Wc now describe in detail the sys- 
tematic effects we have investigated and the mod- 
els we have used for them: 

3.1. Baryonic Infall 

As the baryons within a dark matter halo cool 
and collapse to form a galaxy, the mass distri- 
bution of the halo must respond. This changes 
the depth of the potential well, and thus the 
circular velocity in the halo on galactic scales 
(Blumenthal et al. 1986). Effectively, baryonic in- 
fall remaps a set of dark matter halos to higher cir- 
cular velocity. Although the idea is not new, theo- 
retical investigations of baryonic infall are still im- 
proving (Bullock et al. 2001); we expect that they 
will be significantly advanced by the completion of 
DEEP2. Wc thus have simply adopted a straight- 
forward model from the literature (Gonzalez et al. 
2000). 

That model defines a remapping of the veloc- 
ity function which combines an analytic treatment 
of baryonic infall with results of N-body simula- 
tions. There arc three parameters in this model: 
A, the galaxy spin parameter; Cvir , the concentra- 
tion index for the Navarro, Frenk & White dark 
matter halos assumed by Gonzalez et al.; and 
nid, the fraction of baryonic mass that forms the 
disk. The spin parameter A is well-determined 
by N-body simulations; Gonzalez et al. assume 
A=0.04, which yields results essentially identical 
to integrating over the GIF A distribution (Frenk 
et al. 2000). The concentration index c-uir may 
also be determined from N-body simulations or 
semi-analytically (Bullock et al. 2000); further- 
more, the effects of baryonic infall are essentially 
independent of c„ir for Cyir ^ 5, a range DEEP2 
galaxies at ~ 1 should all fall within. Based 
on N-body simulations, Gonzalez et al. give a re- 
lation between Cyir and Vdrc for 2; = 0, Cyir ~ 
13 X -^Z Voire f^OO km s~^; we take this as our guide, 
with prefactor taken to be 8 rather than 13 (follow- 
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ing the results of Bullock ct al. 2000). This choice 
for our toy model is nearly arbitrary; appropriate 
simulations would be used for baryonic infall cor- 
rections to real data, but we do not require that 
for now. 

However, the final component of the model, m^, 
cannot be determined so reliably from simulations 
or semi-analytic models. We thus adopt it as the 
free parameter in our baryonic infall model, or 
more specifically, mo, the prefactor in the rela- 
tionship between m<j and circular velocity (taken 
to be 0.1 by Gonzalez et al.). This parameter 
has a much stronger impact on the baryonic in- 
fall remapping than A or Cyir, and will be least 
well-known from simulations. The impact of mo 
on baryonic infall-corrected velocity functions is 
shown in Fig. 3. 

3.2. Velocity Errors 

Observations of circular velocities of DEEP2 
galaxies will, of course, be subject to measure- 
ment error. Furthermore, optical spectroscopy 
of emmision-line objects will only determine the 
velocity of ionized gas at some radius within a 
galaxy; it is unlikely that what is measured will be 
identical to the dark halo circular velocity. For in- 
stance, Kobulnicky & Gebhardt (1999) found that 
for galaxies at 2: ~ 0, velocity measurements us- 
ing 0[II] and HI for the same galaxy have a ^20% 
scatter. Because the velocity function follows a 
very steep power law {n{v)dv ~ i;~^) for galaxy- 
scale halos, such a scatter can greatly change its 
shape and therefore potentially affect dN/dz mea- 
surements. For the purposes of this paper we wish 
to find if the amount of scatter may be determined 
from the observed velocity function alone. 

We have adopted a simple toy model for the 
error in velocity measurements: 
ay = a/(25 kms-i)2 + (fvdrc)^- The first term 
arises from the instrumental resolution of DEIMOS, 
which provides a fundamental limitation on the 
precision with which we can measure small ve- 
locities. The second term represents everything 
else that causes scatter between observed veloci- 
ties and Voire (after correction for baryonic infall). 
The results of Kobulnicky et al. suggest that the 
fractional error / ~ 0.1 locally (as the 20% scat- 
ter they found includes the effects of resolution); 
/ ~ 0.2 might be reasonable for galaxies at 2; 1. 
The effect which various values of / have upon the 



velocity function is shown in Fig. 4. 

3.3. Incompleteness 

Two types of incompleteness will affect at- 
tempts to use observable galaxies to count dark 
matter halos halos: incompleteness due to galax- 
ies of a certain luminosity falling beyond the mag- 
nitude limit of the survey, and incompleteness due 
to dark halos which do not contain an observable 
galaxy. The former (which we will hereafter re- 
fer to as "luminosity incompleteness") is easy to 
understand: given a luminosity-linewidth relation 
and its scatter, one can directly calculate the frac- 
tion of galaxies for a given linewidth that will fall 
beyond the magnitude limit of the DEEP2 survey, 
and adjust the velocity function accordingly. On 
the other hand, knowing what fraction of dark ha- 
los with a certain circular velocity host galaxies 
above the surface brightness limit of a survey (or, 
for that matter, contain a galaxy at all) requires 
a full model of galaxy formation. It is not clear 
if this uncertainty, which we will refer to as "halo 
incompleteness" , is an issue for the ^ L^, galax- 
ies DEEP2 will observe or not. Unobserved low- 
surface brightness galaxies could also be a problem 
for the z ^ comparison samples required for a 
dN/dz measurement. However, the local surveys 
now underway will contain 5 — 20 times as many 
galaxies as DEEP2, determining the impact of in- 
completeness from the observed velocity function 
should if anything be more effective than at 2; 1. 

Lacking a model for halo incompleteness, we 
adopt luminosity incompleteness as our toy model. 
Although Bullock et al. (2001) suggest that this 
effect will compromise measurement of dN/dz 
using DEEP2 because the TuUy-Fisher relation 
evolves with redshift, that seems highly unlikely 
to be a problem. With a sample of 10'* galaxies 
that have both magnitudes and linewidths, defin- 
ing the slope, zero point, scatter, and evolution of 
a luminosity-linewidth relation should all be pos- 
sible. The effect of luminosity incompleteness on 
the velocity function is then completely fixed, and 
may be corrected for directly. We demonstrate 
below that in fact we could accurately determine 
luminosity incompleteness (in the absence of halo 
incompleteness) using the velocity function alone. 
However, since we lack sufficient simulations, we 
presume that the impact of halo incompleteness 
on the velocity function would be generally simi- 
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lar, so that we can assess the impact of halo in- 
completeness by evaluating how well we can de- 
termine luminosity incompleteness from the ve- 
locity function alone. This should be gcnerically 
true; it is likely that all high-circular velocity halos 
will contain a galaxy, and one above any surface 
brightness limitations of the survey at that, while 
smaller halos may well not. Thus, the signature 
of the two sorts of incompleteness is similar: full 
completeness for high- velocities halos, rolling over 
to incompleteness at lower velocities. The issues 
involved in halo incompleteness should be much 
better understood by the time the DEEP2 sur- 
vey is completed, as simulations and semianalytic 
techniques should only improve. 

For our toy model of incompleteness, we follow 
the results of Vogt et al. (1997) that the slope 
and scatter of the rest-frame B Tully-Fisher rela- 
tion appear to remain unchanged to z ^ 1. We 
take as a free parameter the zero point of the rela- 
tion, or equivalently W50, the circular velocity (af- 
ter baryonic infall) at which observations to the 
DEEP2 magnitude limit {Iab ~ 23.5) will be 50% 
complete at ^; = 1. We thus obtain a Tully-Fisher 
relation: 

Mb = -7.48 logio + 40.23 - C, (3) 

where C combines the K correction (small for con- 
version of B at z~0 to / at z 1), the Hubble 
constant (as —5 logj^g h) and evolution of the Tully- 
Fisher zero point; the results of Vogt et al. suggest 
C ~ 0.2 ± 0.3 for anh^ 0.75 LCDM model, cor- 
responding to W50 ^ 160 ± 15 km s~^. They found 
that the 0.65 mag dispersion in the Tully-Fisher 
relation at high redshift is consistent with being 
due to a 0.4 mag intrinsic scatter, a 0.2 mag scat- 
ter due to magnitude measurement errors, and a 
0.47 mag scatter due to velocity errors. We apply 
our model incompleteness to the velocity function 
before adding noise, so only the first two elements 
of the scatter are relevant. The effect of varying 
t;50 on the velocity function is shown in Fig. 5. 

3.4. Resulting Errors 

Given these toy models, we have generated 
Monte Carlo realizations of the resulting veloc- 
ity functions as the parameters are varied, as de- 
scribed above. The net effect of these system- 
atics in a standard scenario (mo=0.1, / — 0.2, 
U50 = 160 km s~^) is shown in Fig. 6. We expect 



that one-fourth to one-half of the 50,000 galax- 
ies for which DEEP2 should obtain redshifts will 
provide velocity measurements. We therefore have 
used from 5,000 galaxies (corresponding to the to- 
tal number of objects in the deeper 3HS) to 20,000 
galaxies in our Monte Carlo realizations of the ve- 
locity function. Based upon the simulations, we 
have compiled error budgets for a determination 
of the volume element in three scenarios: a pes- 
simistic one in which there are only 5,000 objects 
with linewidth information and the free param- 
eters of the systematic effects take values which 
result in the greatest errors; a "best bet" scenario 
in which there are 10,000 objects and the free pa- 
rameters have intermediate values; and an opti- 
mistic scenario in which there are 20,000 objects 
and the free parameters take favorable values. We 
present the resulting error budgets for a measure- 
ment of apparent density within the entire DEEP2 
survey volume in Table 1.^ In every case, the dis- 
tributions of the residual errors after correction for 
systematics showed no major non-Gaussianities. 

4. Cosmological Constraints 

The DEEP2 Redshift survey will not simply 
measure dN/dz at one redshift, but instead over 
the entire range 0.7<z<1.5. To determine the 
resulting constraints upon cosmological parame- 
ters, we must take this into account explicitly. 
Because both cosmic variance and the residual er- 
rors from systematic effects after correction are 
essentially Gaussian, we may use a paradigm. 
To simplify, we presume that DEEP2 measures 
dN/dz in eight redshift bins, each covering 0.1 
in z. The observed abundance in each bin will be 
affected by Poisson variance, by the residual er- 
rors from systematic effects after correction, and 
by cosmic variance. 

To determine the Poisson variance in each bin, 
we must know the number of galaxies within 
it. To estimate this, we use a fit to the pho- 
tometric redshift distribution of galaxies in the 
Hubble Deep Fields to the IHS magnitude limit 
(Gwyn 1999; Lanzetta et al. 1999), which yields 

^We have not considered the degree to which systematics 
may be measured or eliminated by study of the local veloc- 
ity function, making our error budget a conservative one. 
Baryonic infall models, in particular, could be further con- 
strained by 2 ~ observations, reducing our dominant sys- 
tematic uncertainty. 
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dN/dz cx (1 + and the total number of 

galaxies appropriate to each error scenario. To 
first order, residual systematic errors will affect all 
redshift bins in the same way. Thus, for these er- 
rors we use the same fractional error in dN/dz for 
all redshift bins and presume it is completely cor- 
related amongst them. Finally, given a cosmolog- 
ical model, we may calculate the cosmic variance 
within each redshift bin as in section 2; we treat it 
as completely independent between bins. In actu- 
ality, the number of objects in adjoining bins will 
be slightly correlated, but this is negligible com- 
pared to the uncorrelated part of the variance (or 
to the residual errors from systematics; the cor- 
related part of the cosmic variance between our 
redshift bins is less than 10% of the total). 

Given these definitions, the covariance matrix 
for our redshift bins is completely determined; we 
list the values used (for an LCDM model) in Table 
2. We may then calculate between any model 
and some nominal, "true" model (e.g. LCDM: 
= 0.3, ng = 0.7, w = -1).3 Observed re- 
sults should be distributed as with two degrees 
of freedom, so the results can be directly trans- 
lated into statistical confidence contours. In Figs. 
7-9 we show the results for all three error budgets. 
Fig. 10 shows the separate effects of cosmic vari- 
ance and residual systematics. For the optimistic, 
best bet, and pessimistic scenarios, the error in a 
measurement of w = —0.7 will be 0.069, 0.086, or 
0.11 if rim is known to ±0.025, as compared to 
0.051 if Poisson errors dominated. The worse the 
errors on fi^, the less relative degradation cosmic 
variance and residual systematics will cause: if the 
error in fi^ were 0.01, the resulting variance would 
be 2.2 — 6.5 X as large as predicted from Poisson 
statistics, while if it were 0.05, precision would bo 
only 15 — 60% worse than the Poisson prediction. 

In Fig. 11, we present the best bet contours 
along with predicted or recent results for a vari- 
ety of cosmological tests. Comparison to Fig. 3 
of NDOO shows that even an optimistic but real- 
istic error budget would yield poorer constraints 
than we previously predicted based on purely sta- 
tistical errors. On the other hand, the measure- 

^We use the extension of to a multivariate distribution 
with covariance: = ~ no)'^V^^(n — no), where n is 
the vector of observations, no is the vector of true values, 
and V is the covariance matrix for no . The elements of V 
are listed in Table 2. 



mcnt should yield much stronger limits on the 
dark energy than any other method available on 
the same timescale, at minimal marginal cost as 
it is a byproduct of a ground-based redshift sur- 
vey that has been designed to address a number 
of issues in galaxy evolution and cosmology. If 
SNAP results are dominated by statistical uncer- 
tainties, they would yield cosmological constraints 
a factor of 2 — 3 x stronger than a dN/dz measure- 
ment based on the DEEP2 galaxy velocity func- 
tion alone. In the long term, the DEEP2 results 
would remain a valuable check on SNAP and other 
methods, however, as they will be subject to very 
different systematic errors than other cosmolog- 
ical tests. DEEP2 will also provide other com- 
plementary measurements of cosmological param- 
eters through a variety of methods (e.g. by study- 
ing the abundance and velocity function of galaxy 
clusters within the survey volume; see Newman et 
al. 2001). 

In summary, we have shown that the uncertain- 
ties of Newman & Davis (2000) were overly opti- 
mistic by a factor of ~ 2, as we ignored the effects 
of cosmic variance and had not yet tested the as- 
sumption that systematic effects could be found 
and corrected for based on the observed velocity 
function. We have demonstrated that this pro- 
cedure can be implemented, though the residual 
errors from the correction procedure somewhat ex- 
ceed the Poisson errors in a dN/dz measurement. 
However, we have found no evidence that the test 
is compromised either by cosmic variance or any of 
the systematic effects considered. We must note 
also that our analysis has been conservative; for 
instance, we have not assumed that observations 
of the velocity function at 2: ~ provides any con- 
straints on systematic errors, nor have we consid- 
ered the fact that comparison to larger-area mul- 
ticolor imaging or redshift surveys could allow us 
to normalize out cosmic variance in the DEEP2 
fields. With sufficient attention to detail, the clas- 
sical dN/dz test can provide much stronger con- 
straints on dark energy than have ever been avail- 
able before. 

We wish to thank Michael Turner, whose chal- 
lenge to produce a realistic error budget inspired 
this work, and Richard Ellis, whose inquiry about 
the importance of cosmic variance caused us to in- 
vestigate the issue. Andrew Jaffe kindly provided 
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plotted here. We also thank Alison Coil, Christian 
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by the National Science Foundation under Grant 
No. AST-0071048. This work was also made pos- 
sible by equipment donated by Sun Microsystems. 
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Fig. 1. — The fractional cosmic variance in a density measurement as a function of the area on the sky- 
surveyed per field (4 fields assumed) in a volume extending from 0.7 < z < 1.5. Curves are plotted for axis 
ratios b/a= 1/8, 1/4, 1/2 and 1. The planned geometry for the DEEP2 survey is indicated with a diamond. 
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Fig. 3. — (top) The PDF for the circular velocity of dark matter halos at ^ = 1 as the disk fraction parameter 
mo is varied from 0.05 (dotted) to 0.2 (dot-dot-dot-dash), (bottom) The ratio of the PDF for varying vahies 
of mo to that in the fiducial model. Also plotted are expected Poisson error bars for that model if the velocity 
function is divided into seven bins. Compare to Figs. 4 and 5; the signatures of the three systematic effects 
differ strongly. 
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Fig. 4. — (top) As Fig. 3, but now the velocity error parameter / is varied from 0.1 (dotted) to 0.4 (dot- 
dot-dot-dash), (bottom) The ratio of the PDF for varying values of / to that in the fiducial model. The 
signature is much larger than the expected errors. 
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Fig. 5. — (top) As Fig. 3, but now the 50% incompleteness velocity V50 is varied from 130 (dotted) to 190 
(dot-dot-dot-dash) km s~^. (bottom) The ratio of the PDF for varying values of V50 to that in the standard 
model. Again, the signature is quite strong. 
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Fig. 6. — The differential circular velocity function for dark matter halos at z = 1, before (dotted) and 
after (solid) the systematic effects described in § 3 are applied. For this plot, the standard values of the free 
parameters (mo = 0.1, / = 0.2, v^o = 160 km s~-^) were used. 
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Fig. 7. — Black curves: 95% confidence constraints in the — plane resulting from our three error 
scenarios in an Qrn — 0-3, = 0.7 model: optimistic (dotted), "best bet" (solid), and pessimistic (dashed). 
Also plotted for comparison is the target 95% statistical uncertainty for the SNAP project (grey, dot-dashed 
curve). 
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Fig. 8. — Black curves: 95% confidence constraints in the ^Im — w plane resulting from our three error 
scenarios in an — 0.3, = 0.7, w) = —1 model: optimistic (dotted), "best bet" (solid), and pessimistic 
(dashed) . Also plotted is the target 95% statistical uncertainty for the SNAP project for comparison (grey, 
dot-dashed curve). 
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Fig. 10. — The best bet contours from the previous three figures (sohd black curves), along with curves 
indicating what the constraints with zero systematic errors (dotted curves) or zero cosmic variance (dashed 
curves). (Top) Constraints in the il„i — plane for an = 0.3, Ha = 0.7 model. (Middle) Constraints 
in the — w plane for a model with = 0.3, w = —1. (Bottom) As the middle panel, but for a model 
with IV = —0.7. In most cases, the effects of both cosmic variance and systematic effects are of comparable 
importance. 
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Fig. 11. — The best bet contours from the previous three figures (solid red curves), along with other current 
or near-future measurements. (Top) Plotted constraints in the ^m — ^A plane include the current 68% SNe la 
results (Perlmutter et al. 1999; large black dashed curve) and BOOMERANG/MAP measurements (Jaffe et 
al. 2000; blue dotted curve), along with (for an = 0.3, Ha = 0.7 model) the 95% SNAP target statistical 
uncertainty (black dot-dash curve). (Middle) Potential measurements in the ilm — w plane for a model with 
rim = 0.3, w = —1. In addition to our "best bet" scenario, the SNAP target statistical uncertainty is shown 
(black dot-dashed curve). (Bottom) As the middle pg^el, but for a model with w = —0.7. 



Table 1 

Error Budgets for Density Measurements from DEEP2 



I Error Scenario: 

Pessimistic Best Bet Optimistic 

Error Source Fractional Error (Parameter Value) 



Counting Statistics 


0.014 


(N= 


5000) 


0.010 


(N= 


10000) 


0.007 


(N= 


20000) 


Cosmic Variance 


0.022 




= 0.723) 


0.018 


(o-8 


= 0.600) 


0.014 


(o-8 


= 0.474) 


Baryonic Infall^ 


0.080 


(mo 


= 0.05) 


0.050 


(too 


= 0.1) 


0.033 


(too 


= 0.2) 


Velocity Errors^ 


0.012 


(/ = 


= 0.4) 


0.002 


(/ = 


: 0.2) 


0.001 


(/ = 


= 0.2) 


Incompleteness ^ 


0.021 


(^^50 


= 175) 


0.014 


(^50 


= 160) 


0.009 




= 145) 



^Residual error when the observed velocity function is used to measure and remove the 
effect 
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Table 2 
CovARiANCE Matrix Elements 



Element^ 


Pessimistic 


Best Bet 


Optimistic 


1,1 


0.00313 nf 


0.00620 nl 


0.01203 nl 


2,2 


0.00307 


0.00611 ni 


0.01196 ni 


3,3 


0.00303 ni 


0.00608 nl 


0.01197 n§ 


4,4 


0.00303 n| 


0.00609 n| 


0.01205 nl 


5,5 


0.00304 ni 


0.00614 nl 


0.01220 ni 


6,6 


0.00308 nl 


0.00623 ni 


0.01241 ni 


7,7 


0.00313 


0.00635 n? 


0.01267 n? 


8,8 


0.00320 n| 


0.00649 n| 


0.01298 n| 




0.00123 rijnj 


0.00305 n^nj 


0.00719 ninj 



""^ These are elements of the covariancc matrix for n^, 
where n^ is the number of objects in the ith redshift bin. 
An LCDM model with = 0.3, = 0.7, w = -1 has 
been used. The first bin extends from z = 0.7 to 0.8, the 
second from 0.8 to 0.9, etc. If fcv is the fractional error 
from cosmic variance in a bin, f count that from counting 
statistics, and fsys that from residual systcmatics (com- 
bined in quadrature), the diagonal elements of the covari- 
ance matrix will be {f^y + /2^„„t + Jn?, while the off- 
diagonal elements will be fsygninj. 
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